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This paper optimises projected land-use structure in 2020 with the goal of increasing terrestrial ecosystem 
carbon storage and simulates its spatial distribution using the CLUE-S model. We found the following: The 
total carbon densities of different land use types were woodland > water area > cultivated land > built-up 
land > grassland > shallows. Under the optimised land-use structure projected for 2020, coastal Jiangsu 
showed the potential to increase carbon storage, and our method was effective even when only considering 
vegetation carbon storage. The total area will increase by reclamation and the original shallows will be 
exploited, which will greatly increase carbon storage. For built-up land, rural land consolidation caused the 
second-largest carbon storage increase, which might contribute the most as the rural population will 
continue to decrease in the future, while the decrease of cultivated land will contribute the most to carbon 
loss. The area near the coastline has the greatest possibility for land-use change and is where land 
management should be especially strengthened. 

Land-use change is a major driver of terrestrial ecosystem carbon storage 1 . Changes of land-use type from 
high- to low- vegetation biomass release carbon into the atmosphere. This not only reduces carbon storage in 
vegetation directly 2 but can also determine the amount of vegetation residue returned to the soil, which is the 
main source of soil carbon storage (SOC) 35 . As high-vegetation biomass, woodland usually has high SOC 
densities compared with other land-use types 6,7 . Changes from woodland to other land uses usually release 
carbon from soil, especially the upper soil 8 . However, there are studies showing cultivated land with higher 
SOC densities than woodland 9 , which may be attributed to the use of organic fertiliser and straw that is returned 
to the soil. Because vegetation and soil carbon densities vary with land- use type, land-use change can affect carbon 
storage, although this may take a long time for soil 10 . With continued urbanisation, increasing areas of cultivated 
land, woodland and wetland have been converted to built-up land, which greatly influences carbon storage in 
terrestrial ecosystems. The impact on carbon balance of land-use and land-cover changes in terrestrial ecosystems 
has become a focus of global change research in recent decades 11 . 

Many scholars have proposed increasing carbon storage by adjusting the land-use structure, which appeared to 
be an effective plan. Lai (2010) 12 and Chuai et al. (2013) 1 optimised land-use structure to increase carbon storage 
at the national and provincial scales in China, respectively, finding that such measures can help China reduce 
carbon emissions effectively. Zhong et al. (2006) 13 optimised land-use readjustment to reduce carbon storage loss 
in Cuiyuan Village in Hubei Province, China. Coastal zones are natural transition regions between terrestrial 
ecosystems and ocean ecosystems. These zones usually have marine characteristics, such as wetlands, shallows 
and salt pans along coastlines. Moreover, 60% of the world's population is concentrated in coastal areas, where 
there is typically a highly developed economy with frequent land-use change 14,15 . The carbon emissions caused by 
such change are greater than in inland areas 1,13 . Some studies have targeted carbon storage variations brought 
about by land-use change in coastal regions, but only in single ecosystems. These studies involved carbon 
emissions caused by deforestation 16 and degradation of wetlands 17,18 . However, no land-use structure optimisa- 
tion based on high carbon storage has been performed for coastal regions. Furthermore, all research examining 
measures to increase carbon storage in other regions has only performed optimisation from the aspect of land-use 
area. Where will the greatest land-use changes be? Which regions will have greater pressure to increase carbon 
storage, and what will be the future spatial distribution of the principal land-use transfers? All these questions 
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Table 1 Carbon d 


ensities of different 1 


and-use types in coastal Jiangsu (t/km 2 ) 












Land use type 






Carbon density 


Cultivated land 


Woodland 


Wetland Water area 


Shallow 


Built-up land 


Vegetation 

Soil 

Total 


543 
8560 
9103 


3360 
8710 
12070 


488 0 
6040 9090 
6528 9090 


107 
5410 
5517 


482 
8267 
8750 



remain unanswered, so the study of land-use structure optimisation 
based on terrestrial ecosystem carbon storage and its spatial simu- 
lation in coastal regions is necessary and meaningful. 

In this study, we use land- use data, soil sample data, and data 
related to vegetation carbon densities and developmental policies 
in the coastal region of Jiangsu Province. Based on these data, this 
work investigates the carbon densities of different land-use types, 
establishes a land-use structure optimisation model to maximize 
terrestrial ecosystem carbon storage, simulates the spatial distri- 
bution of the optimised land use structure in 2020 by CLUE-S model, 
and analyzes the main land transfer types that can affect carbon 
storage and their spatial distributions between 2010-2020. 

Results 

Carbon densities of various land uses and changes. Table 1 shows 
that carbon densities vary widely with land-use type in coastal areas of 
Jiangsu Province. Differences in vegetation carbon densities were 
great, with a highest density of 3,360 t/km 2 and a lowest density of 
zero. The soil carbon density ranged from 5,410 to 9,090 t/km 2 . 
Woodlands have greater total carbon densities relative to other 
land- use types and, in our study, woodland had the highest 
vegetation carbon density and second-highest soil carbon density. 
The total carbon density of cultivated land was similar to that of the 
water area, both of which were much higher than that of wetland; 
neither the vegetation nor soil of wetland had high carbon densities in 
our coastal study area. The carbon densities of shallows were the 
lowest compared with other land- use types, but built-up land had 
high densities, only slightly lower than cultivated land and water area. 

According to the total carbon density values for different land- use 
types, a transfer from cultivated land always resulted in carbon stor- 
age loss, except for the transfer to woodland. Transfer from wood- 
land, with the highest carbon density, to any other land-use type 
caused high carbon loss. This was especially so with the transfer to 
shallow soil, with 6,553 t carbon loss per km 2 , the maximum change 



observed. Wetland loss always increased carbon storage, except when 
a transfer was made to shallow soil. Transfers from water area and 
built-up land to cultivated land and woodland increased carbon 
storage. Transfer from shallow soil, with the lowest carbon density, 
to any other land-use type increased carbon storage. 

Effect of optimised land-use structure on increased carbon storage. 

As indicated in Table 2, the total carbon storage in coastal Jiangsu in 
2010 was 29,759.07 X 10 4 t. Within ten years, variations in carbon 
storage mainly resulted from vegetation changes in new land-use 
areas. Carbon in soil changes slowly and is primarily determined 
by the previous land-use type, so we only calculated the carbon 
storage change due to vegetation. Thus, the regional total carbon 
storage will reach 29,938.75 X 10 4 t in 2020 under the optimised 
land-use structure proposed in this paper. This optimal structure 
will increase carbon storage by 179.68 X 10 4 t compared with the 
2010 structure. A total of 157.26 X 10 4 t of the calculated increased 
storage resulted from reclamation of the new increased area, 
consisting primarily of reclamation as cultivated land (108.40 X 
10 4 t). The remaining 22.41 X 10 4 t of increased carbon storage 
resulted from land-use structure adjustment of the original 2010 
land areas. Without considering the new increased land area, most 
ecological land area will increase, especially the water area, which will 
increase by 893.53 km 2 . However, because we only consider 
vegetation carbon storage and the water area has no vegetation, any 
increases of that area will not increase carbon storage. Wetland and 
woodland will increase by 358.39 and 85.51 km 2 , respectively, but 
because the vegetation carbon density of woodland is much higher 
than that of wetland, the increased woodland area can increase 
carbon storage more. To meet the demand for economic 
development, the decreasing trend in cultivated land area cannot 
readily be stopped and will lead to 21.86 X 10 4 t carbon loss. The 
expansion of built-up land will be controlled effectively, with only 
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Cultivated land 


24325.31 


23922.78 


25919.04 


-402.53 


1593.73 


1996.26 


Woodland 


31 1.30 


396.81 


396.81 


85.51 


85.51 


0.00 


Wetland 


520.47 


878.86 


1340.27 


358.39 


819.80 


461.41 


Water area 


2300.61 


3194.14 
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893.53 


893.53 


0.00 
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1 149.66 


0.00 
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-1 149.66 
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22121.47 
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375.74 
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40.01 
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Water area 
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2091.25 


2091.25 


0.00 


0.00 


0.00 


Shallow 
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-4.70 


18.74 
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4174.72 


4185.07 


4203.82 


10.35 
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7.61 


Total 


29759.07 


29781.49 


29938.75 
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214.76 km 2 increased area, which can increase carbon storage by 
10.35 X 10 4 t in the form of vegetation. 

To clearly analyse the carbon storage change caused by original 
land-use transfer without consideration of the new increased land 
area, we constructed a land- use and carbon storage transfer matrix 
between 2010 and 2020, based on the original study area in 2010. As 
indicated in Table 3, 2,236.13 km 2 land area will be transferred to 
other land use types, representing 7.18% of the entire land area in 
2010. All original shallows will be transferred, with a large area of 
1,149.66 km 2 , accounting for 51% of the entire transferred area. This 
can increase vegetation carbon storage by 29.05 X 10 4 t. The main 
transfers are to wetland, water area, built-up land and cultivated land, 
with vegetation carbon storage changes of 12.33, —3.14, 11.68 and 
8.47 X 10 4 t, respectively. Cultivated land will contribute 786.95 km 2 
land area for transitions. This represents 35% of the entire trans- 
ferred area and will decrease carbon storage by 7.37 X 10 4 t, with 
transitions to built-up land and water areas contributing most to this 
decrease. Transitions to woodland will most inhibit this decrease. 
Because of the consolidation of rural residential land, some built- 
up land has the potential to be converted to cultivated land 
(190.24 km 2 ) and water area (93.92 km 2 ), but this would weakly 
increase carbon storage, because water area has no vegetation. No 
decreases were allowed for woodland, wetland or water area because 
of their high carbon density to maximum carbon storage ratio. 

Spatial distribution of land-use conversion between 2010 and 
2020 and reclamation of new increased area. According to the 
effective CLUE-S model we tested, we produced spatial 
distributions of optimised land-use structure. To discover where 
clear land-use transfer will occur, what the principal transfer types 
will be, and which regions will have greater pressure to increase 
carbon storage, we also produced a land-use transfer map 
comparing 2010 and 2020. Using this and the reclamation plan for 
the new increased area, we constructed Fig. 1. The figure shows that 
coastline regions have a strong potential for land transition, 
exhibiting concentrated distributions for both the original land 
area and the new increased area. Reclamation of the new increased 
area will mainly occur along the southern coastline where there are 
obvious sedimentation effects. Inland regions have dispersed 
distributions, and transitions will mainly occur southwest of 
Yancheng and Nantong and north of Lianyungang, regions where 
cities and towns are mostly located. 

Except for the new increased area, we established six buffer areas 
with distances to the coastline of 0-10 km, 10-20 km, 20-40 km, 
40-60 km, 60-80 km and 80-100 km, to show regional differences 



in land transitions. As indicated in Table 4, the highest percentage of 
transferred land relative to the total buffer area was in buffer area 1. 
This value was 26.9%, much higher than in other buffer areas, with a 
total area of 1,518.26 km 2 , representing 67.32% of the entire trans- 
ferred land area. This indicates that land in the 0-10 km range from 
the coast has the highest probability of transfer. The main transfer 
types include cultivated land to water area; shallow soil to wetland, 
built-up land, water area and cultivated land; and cultivated land to 
built-up land. In buffer 2, the percentage decreased sharply to 4.08%, 
with main transfers of cultivated land to water area, built-up land and 
wetland. Buffer 3 showed the lowest percentage and smallest absolute 
transfer area at 98.59 km 2 . Buffers 4, 5 and 6 showed increasing 
percentages from 1.73% to 9.13%, with absolute areas from 137.53 
to 106.26 km 2 . The much higher value of 9.13% in buffer 6 is attrib- 
utable to widespread water area there; the nearby land is under great 
pressure to be replaced by water area. 

Without consideration of soil organic carbon, we produced a map 
of vegetation carbon density changes in Fig. 2 to show the spatial 
distribution of carbon sinks or sources generated by land-use change. 
This shows that the reclamation of the new increased area will 
strongly enhance carbon storage because the original sandy beach 
will be covered by vegetation, especially in the case of reclamation to 
cultivated land. For the original land area transitions, most transi- 
tions near the coastline have great potential to increase carbon stor- 
age, especially in the central and south regions. Part of the northern 
coastline is under pressure to lose carbon storage capacity, because of 
water area expansion. The transition of cultivated land to built-up 
land in those areas will also lead to carbon loss. For inland districts, 
most of the transferred area will enhance carbon storage. The regions 
of carbon loss are mainly located in western Yancheng. This is prim- 
arily caused by the transition of cultivated land to water area. 

Discussion 

Land-use change caused by human activities is a main driver of 
modifications to terrestrial ecosystem carbon storage 1 . Coastal 
regions usually have highly developed economies, which drive fre- 
quent land-use changes 14,15 . Using land- use type maps from 2005 and 
2010, vegetation-type maps, soil sample data, and other empirical 
data related to carbon densities for different terrestrial ecosystems, 
and based on a linear programming model and the CLUE-S model, 
we examined carbon densities of both vegetation and soil for various 
land-use types. Based on these carbon densities, we optimised land- 
use structure in 2020 based on maximum carbon storage in terrestrial 
ecosystems and simulated its spatial distribution in coastal Jiangsu. 
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Cultivated land 


23538.36 


71.43 


34.40 


480.16 


0.00 


200.96 


24325.31 


Woodland 


0.00 


31 1.30 


0.00 


0.00 


0.00 
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31 1.30 


Wetland 


0.00 


0.00 
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0.00 


0.00 


0.48 


520.47 


Water area 


0.00 


0.96 
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0.00 


0.00 
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Shallow 


1 94. 1 8 


0.00 


323.67 


320.41 


0.00 


31 1.40 


1 149.66 


Built-up land 
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0.80 
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396.81 
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Cultivated land 


0.00 
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Woodland 
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Wetland 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


Water area 


0.00 


0.32 


0.00 


0.00 


0.00 


0.00 


0.32 


Shallow 
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0.00 
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0.00 


1 1.68 


29.05 


Built-up land 
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0.00 
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Figure 1 | Spatial distribution of land-use transfer between 2010 and 2020. Map created using ArcMap 9.3 software. (Environmental Systems Resource 
Institute (ESRI), Redlands, CA, USA). 



Between 2010 and 2020, we used only vegetation carbon storage to 
analyse the effects of the optimised land-use structure on increasing 
carbon storage. Because carbon in soil changes slowly, possibly 
occurring over periods longer than ten years 10 , carbon storage for 
new land-use types is mostly determined by previous types. We also 
assumed that no land-use types gain or lose carbon over time. The 
only changes in carbon storage considered were those resulting from 
the conversion of one land-use type to another. Our study area is 
mainly covered by crops and wetland vegetation. Crop vegetation 
covers most of the area, and it has relatively stable growth because of 
agricultural management. Wetland vegetation is mainly located 
along the coast, and large areas have been effectively protected by 
ecological protection zones 19 . Thus, vegetation growth status was not 
considered in our study, 

Woodland had the highest carbon density in all regions, which can 
be explained by the fact that it always has high levels of biomass 



sequestered as vegetation 12 . Moreover, woodlands create more resid- 
ual vegetation, which can also enhance carbon storage in soil 4 . 
Although lacking vegetation, water areas also showed high SOC 
densities. This is likely attributed to water transport mechanisms that 
can collect organic carbon from other areas and accumulate it in silt 
under the water 20 , as well as to long-term flooding conditions that 
reduce decomposition rates 7 . The carbon density of cultivated land 
was only slightly less than that of water area, which can be explained 
by the fact that cultivated land has a greater vegetation biomass than 
all other land-use types except woodland. For soils, the effects of 
returned crop residues and straw, long-term use of organic fertilisers, 
and flooding of paddy land in our study area can also enhance carbon 
storage 21,22 . Built-up land has a vegetation carbon density similar to 
wetlands. This is due to the high vegetation coverage rate in coastal 
Jiangsu. The high soil carbon density agrees with the results of earlier 
studies 23,24 and can be explained by the multitude of soil carbon 
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t-up land to cultivated land, cultivated land to water area and woodland 
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Figure 2 | Spatial distribution of carbon density changes between 2010 
and 2020 (t/km 2 ). Map created using ArcMap 9.3. 

sources in towns, including vegetation such as leaves, twigs and 
weeds, substantial amounts of household garbage, and organic waste 
produced by urban industries 7 . Additionally, soil in urban areas is 
generally compacted and sealed, which can prevent the release of soil 
carbon into the atmosphere, thereby reducing carbon decomposition 
and increasing its accumulation 7 . The carbon density of wetlands was 
relatively low. This contradicts earlier studies that showed higher 
carbon densities for wetlands compared with other terrestrial eco- 
systems, especially for soil 25,26 . This can be explained by the fact that 
wetlands in our area are mainly composed of halophytic vegetation 
and are located in coastal solonchak soil areas with high sand con- 
tent. These soils do not favour the concentration of carbon in the 
soil 27 . Shallows have little vegetation coverage and are a weak carbon 
source, and the sand content in solonchak soil may be even higher 
than that in wetland areas, so shallows had the lowest carbon densi- 
ties overall. 

The reclamation of new increased area can greatly benefit carbon 
storage. Without considering the new increased area and the SOC of 
the entire area, we can increase carbon storage by 22.41 X 10 4 t in 
coastal Jiangsu by adjusting the land-use structure. This can offset 
1.88% and 1.05% of the carbon emissions resulting from increased 
regional energy consumption in 2005 and 20 10 24 , respectively. If we 
fail to optimise and permit the land-use structure to change as it did 
in 2005-2010, many more ecological land-use types, such as culti- 
vated land, woodland and wetland, will decrease sharply. Moreover, 
the expansion of built-up land will occur more rapidly. The much 
greater decrease of ecological land will lead to tremendous carbon 
loss from both vegetation and soil. The rapid expansion of built-up 
land will also lead to more anthropogenic carbon emissions. Because 
the principal activities of energy consumption are tied to certain 
types of such land, there are usually high carbon emissions from 
energy consumption. Putting limits on this land type can therefore 
obviously decrease carbon emissions from energy consumption 12,28 . 
Ecological lands can also provide benefits other than carbon storage, 
such as water conservation and flood control, erosion control, fuel, 
food, biodiversity and others 29 . Our study did not consider carbon 
changes due to soil. If there is sufficient time for the transferred land 



to mature into another ecological land system, the optimised land 
structure will cause another 314.18 X 10 4 t SOC increase, without 
consideration of the new increased area as a long-term effect. Some 
transition types in our optimised land-use structure that can lead to 
carbon loss may yield opposite effects, such as the transition of cul- 
tivated land to water area. Although this process will not be com- 
pletely finished during the period 2010-2020, there will be real SOC 
changes to some extent during that decade, especially in the surface 
soil layer. Consequently, the actual increased carbon storage may be 
much greater than 22.41 X 10 4 t from vegetation only. Furthermore, 
many land-use types currently undergoing natural succession can 
continue to accumulate carbon for a very long time 30 , because our 
analysis was based on the assumption that none of the land-use types 
gain or lose carbon over time. If vegetation is well protected, carbon 
storage in 2020 will increase because of the carbon sink effect of 
vegetation. Therefore, the optimised land use structure is effective, 
especially over the long-term. 

The spatial distribution simulation indicated that the area near the 
coast had the greatest probability of land- use change. This is deter- 
mined by the development policy for the coastal zone proposed by 
the provincial government. Thus, we should pay more attention to 
the near-coastal area. In addition to land-transfer type control, we 
can also adopt other means to increase carbon concentrations in this 
important region, including land management, tillage practices such 
as combined use of inorganic fertiliser and organic manure on culti- 
vated land, and the application of a suitable crop rotation system 31 . 
The reclamation of sand beaches and shallows can produce most of 
the increased carbon storage, but we should pay attention to the 
protection of the ecological environment during the process of 
exploitation 32 . Because the rural population is declining, the demand 
for rural residences will shrink. Rural residential land will most likely 
be converted to cultivated land or other land types, and this repre- 
sents a high potential for increasing carbon storage. Thus, we should 
strengthen rural land consolidation, which may contribute the most 
to increasing carbon storage in the future. 

In summary, our study may be valuable to land managers and local 
governments, for we examined carbon storage capacities of different 
land use types. Our method can effectively increase carbon storage in 
coastal Jiangsu and can guarantee land-use demands for coastal eco- 
nomic and social development. Our spatial simulation analysis cla- 
rified the spatial distribution of land- transfer types, which can make 
spatial land management more targeted and effective. Therefore, this 
study can be highly beneficial to the ecological and sustainable 
development of coastal Jiangsu. 

Methods 

Study region. The coastal region of Jiangsu faces the Yellow Sea to the east and the 
Yangtze River to the south. The study area encompasses the cities of Lianyungang, 
Yancheng and Natong (Fig. 3), with an area of 3.3 X 10 4 km 2 , or 35% of Jiangsu's total 
land area. The area lies between 31 41'N-35°07'N and 118°24'E-121°55'E. The total 
length of the coastline is 954 km. In recent years, the economy has developed rapidly, 
with frequent land use changes 14 . 

Data sources. The 100 m X 100 m grid land-use map was produced from Landsat 
Thematic Mapper data. The vegetation-type map was compiled using data from the 
2000s and was able to effectively describe the most recent vegetation distribution in 
the study area. There were 7,797 soil samples obtained from a multi-purpose regional 
geochemical survey in Jiangsu Province from 2003 to 2007. Forest data were acquired 
from The Fifth Forest Resource Inventory in the province. Crop yield data were 
obtained from statistical yearbooks of the three coastal cities. Other empirical data 
were obtained from related research studies. 

Spatial driving factor maps included population density, GDP density, elevation, 
and distances to main roads, cities, villages, ports, coastline, and main water lines. 
Population and GDP data were obtained from statistical yearbooks of the three 
coastal cities. A digital elevation model (DEM) was created based on Shuttle Radar 
Topography Mission data provided by the International Scientific Data Service 
Platform. 

All figures (1-3) in our paper were generated by geographic information system 
mapping techniques. 
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Figure 3 | Location of study area. Map created using ArcMap 9.3. 



SOC densities for different land uses. For each soil sample, the SOC density was 
calculated by the following equation: 



SOCD= 1000 x TOC xpxh, 



(1) 



where SOCD represents SOC density, TOC is the soil organic carbon content (%), p is 
the bulk density, and h is the soil depth. 

Geographic information system mapping techniques were used to produce soil 
sample distribution maps according to latitude and longitude. Kriging methods based 
on a spherical model were then used, and the sample point data were converted to 
polygon data covering the entire study area. Based on the land-use map of 2005, zonal 
statistics in ArcGIS 9.3 software (Environmental Systems Resource Institute (ESRI), 
Redlands, CA, USA) were used to compare soil carbon densities under various land 
uses. 

Vegetation carbon densities for different land uses. 

1) Woodland. This was calculated from biomass density and carbon content. 
According to related previous studies, the carbon content of trees is usually 
taken as 0.5 33 , and the calculations are shown below. 



i x 0.5 



= «+fcV», 



(2) 



where W is the average vegetation carbon density of woodland; Wp S« and 
are the vegetation carbon density, woodland area and vegetation biomass 
densities for tree type i and tree-age j, respectively; V« is tree volume density; 
and a b are constants obtained from the study of Xu et al. (2007) 33 , who 
established a linear regression equation between forest biomass and its volume, 
based on 2304 forest sample plots in China. 

Cultivated land. Major crops in coastal Jiangsu include rice, wheat, corn, beans, 
potato, cotton, peanuts, rapeseed, sugar and vegetables. Here we calculate the 
crop vegetation carbon density according to cropland yield as shown below: 



C = 



Q = 



Pi x (1 - Wt) x Y } 
Hi x Si ' 



(3) 



where C is the average vegetation carbon density of cultivated land; Q is the 
vegetation carbon density of crop type i; 5, is the area of crop type i; and P ; , W, 
Y h and ff, represent the crop yield, water content, ability to absorb carbon, and 
the economic coefficient of crop type /, respectively. S, and P, were obtained 
from statistical yearbooks of the three coastal cities, and values for W, Y,- and H, 
were acquired from the results of related research by Zhao (201 1 ) 34 , who carried 
out a similar study in Nanjing, Jiangsu. 

Other land. For wetland, biomass data for various types of vegetation are 
quoted from the research of Zong et al. (1992) 35 , who studied the biomass of 
saltwater vegetation and sand vegetation on the coast of Jiangsu Province. 
Shallows can be divided into Artemisia halodendron-dominated beaches, 
marshes dominated by Spartina and bare beaches. Therefore, according to 



the vegetation type and area, the vegetation carbon density of shallows can 
be calculated. The vegetation carbon density of trees and grass scattered among 
built-up land was determined according to the vegetation coverage rate and the 
mean values of tree and grass carbon density. We defined the vegetation carbon 
density of water area as zero, because it has almost no vegetation coverage. 



Land-use structure optimisation. A linear programming model was selected to 
optimise land-use structure. This model includes the establishment of a target 
function (Equation 4) and constraint conditions. 



T= y2Ai(Vi-\-Si)= max i = l 5 2,3,...,n, 



(4) 



where T is the regional total carbon storage of the terrestrial ecosystem; A { is the area 
of land-use type i; and Vj and S, are the carbon densities of the vegetation and soil of 
land-use type i, respectively. 

We established constraint conditions for six variables as follows: cultivated land X l5 
woodland X 2 , wetland (with herbaceous vegetation) X 3 , water area X 4 , shallows X 5 , 
and built-up land X 6 . Here we set 2010 as the initial year and 2020 as the targeted 
prediction year. 

About a third of the coastline in the study area suffered erosion 36 , mainly within 
Lianyungang. The remaining area suffered deposition effects. This effect was much 
more obvious than erosion, and it expanded the total land area annually. The erosion 
approached the coastal levee. This will weaken or stop the erosion effect 37 , and the 
total study area will continue to increase in the future. In contrast to our previous 
study 24 , here we considered the land reclamation plans of local government. The 
affected land is principally along the coastline and offshore sand beaches, which will 
mainly be treated by engineering efforts. Thus, we predict that the total area will reach 
36,935.83 km 2 by 2020: 



^X, = 36,935.83 X,>0 



(5) 



To guarantee local food self-sufficiency, 21,925.19 km 2 of cultivated land will be 
needed, which we consider a low value. According to the land consolidation potential 
of rural residential and land reclamation plans, we predict that cultivated land can 
reach 25,919.75 km 2 by 2020, which we consider a high value. 



21,925. 19<Xi < 25,919.75 



(6) 



Woodland should be strictly protected, such that its area in 2020 should at least be 
higher than in 2010 (311.3 km 2 ). According to the 12th Five-Year Plan of forestry 
development planning and the demands of ecological construction, the area of 
woodland will increase by 85.51 km 2 in 2020 over that in 2010 with effective imple- 
mentation. The constraint conditions are established as: 



311.3<X 2 <396.81 



(7) 



Considering Yancheng coastal agricultural development planning, we predicted in 
our previous study 26 that wetland will reach at least 687.14 km 2 by 2020, because of 
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planted grass. According to the land reclamation plan, part of the existing shallows 
and new increased sand beaches will be reclaimed as wetland, which would increase 
the area by 653.13 km 2 . Wetland should be strictly protected because of its significant 
carbon sink effect, so we establish the constraint conditions as follows: 

X 3 > 1,340.27 (8) 

According to water conservancy development planning of the three coastal cities, at 
least 94.3 km 2 of new area will be needed compared with 2010, which we consider a 
low value. The increased water area is mainly converted from wetland, and the total 
water area will reach 2,896.35 km 2 according to the rate of increase between 2005 and 
2010. We consider this a high value, because wetland will be protected strictly as 
stated above. Moreover, according to the reclamation plan, another 297.19 km 2 of 
existing shallows will be converted to water area. Therefore, we should add this to the 
low and high water-area values above: 

2,692.69<X 4 <3,194.14 (9) 

According to the reclamation plan, original shallow land in 2010 will be entirely 
converted to other land-use types. However, another 710.81 km 2 of new increased 
shallows will be present, owing to natural and artificial boosting of the deposition 
effect: 

X 5 =710.81 (10) 

According to population, per capita living space and land demands for traffic, at least 
4,855.88 km 2 of built-up land will be needed in 2020, which we consider a low value. A 
high value is that used in land-use planning in Jiangsu Province {5,274.49 km 2 ), 
which mainly considers the demands of economic development. According to the 
land reclamation plan, another 518.88 km 2 of built-up land will be needed, mainly for 
harbour construction, which will be reclaimed from existing shallows and the new 
increased sand beaches. Thus, we must add 518.88 km 2 to the high and low values: 

5,374.76 <X 6 < 5,793.37 (11) 



Spatial distribution simulation of optimised land-use structure. We selected the 
CLUE-S model to perform the spatial distribution simulation of the optimised land- 
use structure. This simulation typically considers the following: (1) Land-use 
requirements (demand). For the various land-use types, these are calculated with the 
linear programming model described above. (2) Spatial restrictions. We chose major 
natural reserves as restrictions. (3) Land-use type conversion settings. According to 
the land-use transfer situation and model simulation accuracy, values of conversion 
elasticity for different land-use types were set as follows: cultivated land 0.5, woodland 

0. 6. wetland 0.3, water area 0.4, shallows 0.5, and built-up land 0.8. We set 2005 as the 
initial year and 2010 as the target year to test the accuracy of the CLUE-S model. The 
model passed the ROC and Kappa index tests well, and it was able to simulate land- 
use spatial distributions in 2020. 
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